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Abstract 

One of the main challenges facing wireless sensor networks (WSNs) is the limited power resources 
available at small sensor nodes. It is therefore desired to reduce the power consumption of sensors 
while keeping the distortion between the source information and its estimate at the fusion centre (FC) 
below a specific threshold. In this paper, given the channel state information at the FC, we propose a 
subset selection algorithm of sensor nodes to reduce the average transmission power of the WSN. We 
assume the channels between the source and the sensors to be correlated fading channels, modeled by 
the Gilbert-Elliott model. We show that when these channels are known at the FC, a subset of sensors 
can be selected by the FC such that the received observations from this subset is sufficient to estimate 
the source information at the FC while maintaining the distortion between source information and its 
estimate below a specific threshold. Through analyses, we find the probability distribution of the size 
of this subset and provide results to evaluate the power efficiency of our proposed algorithm. 

Index Terms 

Wireless Sensor Networks, Correlated Fading Channels, Gilbert-Elliott model 

I. Introduction 

Wireless sensor networks (WSNs) [HI, O, [O, [|4| are receiving increasing attention due 
their numerous current and foreseen applications in several fields including field trials and 
performance monitoring of solar panels [El, target detection through digital cameras [[H, and 
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petrochemical industry fields 01 . One of the main challenges of WSNs is to overcome their 
energy constraint problem. That is, sensors are powered by batteries with limited energy budgets. 
Due to deployment of sensor nodes in inaccessible or hostile environments, recharging these 
batteries is often not an option. Also, the network is expected to have a lifetime in the order 
of several months, or even years |[8l. Therefore, design of power-efficient WSNs is of extreme 
importance. Numerous works are dedicated to the topic of energy conservation in WSNs. A 
comprehensive review of these works is given in |[8l. 

A typical sensor node in a WSN consists of three main subsystems namely; sensing, pro¬ 
cessing, and wireless communication subsystems. These subsystems are responsible for data 
acquisition, local data processing, and data transmission, respectively flU. In addition, a power 
source is included in the sensor node with a limited power budget. In llH, it is shown using 
experimental measurements that in most cases data transmission is the most energy-consuming 
unit of the sensor node. In a similar conclusion, it is estimated in lUOll that transmitting one bit 
by the data communication unit requires an energy equivalent to performing about a thousand 
operations in the data processing unit. It is worth mentioning that in some applications, the 
sensing subsystem might consume more power than the data communication subsystem (see ([H 
for details) but in typical applications of WSNs, the highest portion of power is consumed by 
the data communication subsystem. Therefore, it is highly desired to develop protocols to reduce 
the transmission power of sensor nodes and hence extend their lifetime. 

In this paper, we consider a WSN where the source-sensor channels are correlated fading 
channels modeled as Gilbert-Elliott channels ifTTI . ifT^ . The Gilbert-Elliott model is a first- 
order Markov model for a correlated fading channel quantized to binary levels of Good and Bad 
states by setting a proper Signal-to-Noise Ratio (SNR) threshold. The channel in its Good and 
Bad states is modeled by binary-symmetric channels (BSCs) with crossover probabilities of pc 
and pb, respectively. The Gilbert-Elliott model is the simplest possible finite-state Markov (ESM) 
model for correlated fading channels. The problem of modeling a correlated fading channel by a 
ESM process is considered in numerous works. An excellent review of works on ESM modeling 
of fading channels is provided in IIT3l where the relations between real-valued fading channel 
parameters and the ESM channel parameters are also considered. 

Eet a binary source block consisting of M bits be transmitted to N sensors via independent 
Gilbert-Eliott channels. Eor each source-sensor channel, the channel states during this transmis¬ 
sion can be expressed as an M bit binary sequence where we let a bit 1 represent a Good state 
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and a bit 0 represent a Bad state. We eall this M-bit sequenee as the ehannel-state information 
Sequenee (CSI sequence). For slowly varying fading channels the CSI sequence consists of a 
few runs and is efficiently compressed using a run-length code (See Fig. [3] and Table U)- 

Our main contribution in this paper is to propose and analyze a two-phase transmission scheme 
as follows. At the first phase, each sensor compresses its respective source-sensor CSI using a 
run-length code and transmits it to the FC. Based on the received CSI from all nodes, the FC will 
know the location of Good bits, i.e. the bits that are received in a Good channel state. The FC 
then finds the smallest subset of sensors such that for each source bit, at least one of the sensors 
in the subset has a Good observation of that source bit. In other words, this is the subset with 
minimum number of sensors, such that for each source bit at least one of the sensors received 
this bit in Good channel state. Then, the FC sends a feedback signal to request transmission 
from this subset. Therefore, at the second phase, only a subset of sensors transmit to the FC, 
resulting in reduction in the average transmission power. 

The motivation behind our proposed algorithm is as follows. According to the Gilbert-Elliott 
model dm, [IT2l . we have po < Pb, be- Good bits are more reliable than Bad bits. Therefore, 
we are in fact attempting to find the minimum number of sensors such that if these sensors 
transmit to the FC and the rest of sensors remain silent, the FC still receives one (or more than 
one) reliable copy of each source bit and consequently is able to reliably reconstruct the source 
information. To examine this idea more precisely, let the WSN have an end-to-end distortion 
requirement of D < D, where D is the expected value of the normalized Hamming distortion 
(the Bit Error Rate) between the source and its estimate at the EC; and Z) is a fixed distortion 
threshold. If a (minimum-sized) subset of sensors exists such that each source bit is received 
through a Good channel by at least one of the sensors in the subset, then the EC will be able 
to reconstruct the source with a distortion less than or equal pc- Let u be the probability of 
existence of such subset. Then we could bound the end-to-end distortion of the WSN as D < Du 
where Du = u x pa + {1 — u) x f where we used the fact that in worst case, the distortion 
is bounded by |. In Section IV-Bi we show that for a WSN with sufficiently large number of 
sensors, the value of z/ is arbitrarily close to 1 and therefore, limAr^oo Du = Pg- The value of pc 
could be expressed as po = Pb ('^) / (-^1-^ > dX where A* is the SNR threshold applied for 
quantizing the fading channel, Pb (A) is the bit error probability for SNR of A, and / (A|A > A*) 
is the conditional probability distribution function of the SNR. Assuming a binary-phase shift¬ 
keying (BPSK) modulation and an additive white Gaussian noise with two-sided power spectral 
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density of ^ at the receiver, we have Pb (A) = Q , where Q (.) represents the Q-function. 

From the above results, we could bound pc sls pa < Q (\/2At) where for obtaining this upper 
bound we used the fact that Pb (A) is a decreasing function of A and has its maximum value at Xt. 
In conclusion, for WSNs with sufficiently large number of sensors, the distortion upper bound 
Du is always less than or equal to Q (y/2Xt). Therefore, if Xt is such that Q [y/2Xt) < D, we 
could conclude that our subset selection algorithm satisfies the distortion requirement of D < D, 
while reducing the average transmission power of the sensor nodes. In this paper we assume 
that the condition pc ^ D holds, and proceed with presenting our subset selection algorithm. 

The rest of this paper is organized as follows. In Section HI] we present our system model used 
in the paper. In Section III, we present our proposed two-phase algorithm with some examples. 
In Section jlVl we analytically derive the probability distribution of the size of the minimum-si/e 
subset, as a function of network size, channel parameters, and the source sequence length (the 
size of this subset is a random variable that depends on the CSI realizations). We also consider 
the computational complexity of our analytical solution and provide suggestions to reduce this 
complexity in Section V. In Section jVll we provide numerical results to evaluate the efficiency 
of our scheme in terms of power conservation. Finally, Section IVIII concludes the paper. 

II. System Model 

We consider a data gathering WSN illustrated in Fig. [U where an M-bit binary source is 
sensed by N sensors via Gilbert-Elliott channels and then transmitted to the FC via noiseless 
channels. To justify the assumption of noiseless sensor-FC channels, we note that according 
to IEEE 802.15.4 standard, it is recommended that the network combines cyclic redundancy 
check (CRC) codes with automatic-repeat request (ARQ) and continues re-transmission for a 
pre-determined number of times [fT4ll . Therefore, assuming genie CRCs, a sensor’s data is either 
eventually delivered to the EC error-free, or not delivered to the EC at all. We assume that the 
N sensors of Eig. [U are the sensors that succeeded to deliver their data to the EC before the 
maximum allowed number of re-transmissions is reached. Also, note that several researchers 
suggested including a forward-error correction (EEC) scheme at sensor nodes to reduce the error 
probability of the sensor-EC link (e.g., ifTSll . ifT^ . ifTTll . ifT^ and references therein). This will 
reduce the expected number of re-transmission requests. 

The state diagram of the Gilbert-Elliott channel is shown in Eig. [2] The channel is modeled by 
a Good and a Bad states and at each state the channel acts as a BSC with transition probabilities 
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of pg and pb > Pg {.Pg,Pb < 0.5), respectively. The transition probabilities from the Good 
state to the Bad state and from the Bad state to the Good state are represented by parameters 
e and p respectively. As mentioned in Section HI the Gilbert-Elliott channel can be considered 
as a quantized version of a correlated fading channel. Figure [3] shows an example of channel 
realizations for a network with N = 6 sensors and M = 256 source bits. The dark areas 
show the Good state and the white areas show the Bad state. To obtain these realizations, we 
generated 6 realizations of correlated Rayleigh fading channels using Jakes model [|T9ll . Then, 
we applied a quantization threshold of at = 1 on the fading amplitude, a. The fading channels 
have a normalized fading rate of fdTg = 2 x 10“^ where fd is the Doppler frequency and Tg is 
the symbol period. Through Monte-Carlo simulations, we estimated the resulting Gilbert-Elliott 
channel state transition probabilities as e = 0.0075 and p = 0.0041, respectively. It is observed 
from Fig. [3] that the CSI consists of a few runs and therefore, could be efficiently compressed 
by a run-length code. In Table H] we show the expected value of the compression rate for the 
run-length coding scheme, for slowly varying fading channels with different normalized fading 
rates. 


III. Proposed Two-Phase Transmission Algorithm 

Assume that we wish to re-construct the binary source at the EC with a normalized Hamming 
distortion less than or equal to a threshold, D. Also assume that Pg ^ D and define a coverage 
event as follows: 

Definition: A source bit is covered by a subset of sensors if it is sensed via a Good channel 
by at least one of the sensors in the subset. An M-bit source sequence is covered by a subset 
of sensors if all of its bits are covered by the subset. 

For example, in Fig. [3] the source sequence is covered by the subset consisting of the first, 
second, and fourth sensors. Given the above definition, our proposed transmission scheme is a 
two-phase scheme as follows, (i) At the first phase, the sensors transmit their compressed CSI to 
the FC and then wait for a feedback signal from the FC to proceed, (ii) The FC de-compresses 
the received CSI and selects the smallest subset of sensors that cover the source sequence. The 
implementation of the selection algorithm at the FC is shown in Fig. IH As shown in Fig. IH if no 
subset is covering the source sequence, the FC requests transmission from all N sensors, in order 
to collect all available information for reconstructing the source information. After selecting this 
minimum size subset, by transmitting a limited feedback (e.g. an 77-bit string where the selected 
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sensors are marked by 1 and the non-seleeted sensors are marked by 0) the FC informs the 
sensors of whieh subset is seleeted, and only that subset of sensors transmit their observations to 
the FC. It is clear that receiving observations from this subset is sufficient to recover the source 
information with a distortion less than or equal to pc- If Pg is less than or equal to the tolerable 
distortion threshold of the network, which we represent by D, then the received transmissions 
from the selected subset is sufficient to satisfy the distortion requirement of the system. Also, 
by applying this subset selection method, only a portion of sensors transmit at each time and 
therefore, the average transmission power of sensors reduces. 

Let us refer to the size of the selected subset by K. Obviously, K is a. random variable that 
depends on the CSI realizations and takes values from 1 to N. The expected value of K is an 
important indicator in our proposed scheme. The ratio of this expected value to the total number 
of sensors, N, represents the average ratio of sensors transmitting to the FC. If this ratio becomes 
smaller, the average transmission power is reduced. 

To quantify the power efficiency of our proposed two-phase scheme, we consider the total 
number of transmitted bits by sensors as an indicator of the consumed power, and compare this 
parameter with a conventional one-phase scheme where all sensors transmit all their observed 
bits to the FC and no CSI is transmitted. Let us denote the total number of transmitted bits 
of the conventional scheme and our scheme by Bi and B 2 , respectively. Obviously we have 
Bi = M X N. Also, it is easy to observe that B 2 is a random variable and if the expected 
value of the compression rate of the run-length coding scheme is represented by p then, we 
have E [S 2 ] = M x {p + E [K]). Now, if we define an efficiency factor p as the ratio of Bi and 
E [i? 2 ], we have: 

N 

P = —-pr^. ( 1 ) 

p + E[K\ 

If p is greater than one, then our proposed scheme consumes less power compared to the 
conventional scheme. In Section |Vll we evaluate p for Gilbert-Elliott channels with different 
parameters, as well as for different number of sensors and source sequence lengths, M. Our 
results show that in many cases, p is considerably larger than one. 

IV. Probability Distribution of the Selected Subset Size 

As mentioned in Section III, the size of the selected subset is a random variable that depends 
on CSI realizations. Refer to this subset size by K and let fx (.) and Ex (.) be the probability 
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Fig. 1. System model of the data gathering Wireless Sensor Network. 



Fig. 2. Gilber-Elliott channel model 

mass function and the cumulative mass function of K, respectively. Obviously, Fx (k) is the 
probability that there exists a subset of k sensors to cover the source sequence (if a subset of 
smaller size covers the sequence, we could add arbitrarily selected sensors to this subset to 
make its size equal to k). Also, the probability mass value, fx (k) = Fx (k) — Fx {k — 1) is the 
probability that k is the smallest size of a subset that covers the source sequence. In the sequel, 
we derive an analytical expression for Fx{k). 

We assume N independent Gilbert-Elliott channels between source and sensors. Let {nm ^n) 
represent the state transition probabilities for the channel from the source to the sensor number n. 
Assume the transmission of bit number m for a fixed m. Let us denote the source-sensor channel 
states at time interval m by Cm = {Cm (1), Cm {N)) where Cm {n) = 1 if the channel from 
the source to the nth sensor is in Good state, and Cm {n) = 0 if the channel from the source to 






Fig. 3. Realization of Gilbert-Elliott Channel State Information for 6 sensors. The Gilbert-Elliott channels are results of quantizing 
correlated Rayleigh Ending channels with normalized fading rates of fdTs = 0.002. A fading amplitude of at = 1 is used as 
the threshold for quantization. Dark areas show Good channel states (i.e. amplitudes above the threshold) and white areas show 
Bad channel states. 



Eig. 4. The sensor selection algorithm at the Fusion Centre. 


























the nth sensor is in Bad state. Let S be an ordered subset of (1,2, ...,N) with eardinality [S'! 
sueh that = (5(1), 5(2),..., 5(|5|)) and 5(1) < 5(2) < ... < 5(|5|). Define: 


7„.(S) 


1; n’:.=,(E!fi.C'„.(S(n)))>0 

0; otherwise. 


( 2 ) 


In ([2]), 7 m(5) = 1 if at every bit interval m' = 1 : m, the ehannel state from the souree to 
at least one of the sensors in set 5 is in Good state, i.e. all m bits are eovered by the set 5. 
Note that if at some bit interval m' all these ehannel states are Bad, then for that bit interval 
Xlfii Cm' (5 (n)) = 0 whieh results in 7m(5) = 0. 

Using the above definition, Fk (k) is equal to the probability that there exists at least one set 
5 with |5| < /c sueh that 7 m( 5) = 1 (M is the total number of transmitted souree bits). For 
calculating this probability, it is sufficient to calculate the probability that there exists a set with 
|5| = /c and 7 m( 5) = 1 (as mentioned above, if a set with cardinality less than k covers all bits 
up to bit M, we could add a proper number of arbitrarily chosen sensors to make the cardinality 
of this set k and the extended set still covers all bits up to bit M). 


Define 


N 

k 


Now we can write: 


. There exist sets 5 with |5| = /c which we refer to as si, S 2 , 


Fxik) = P((7m(si) = 1) II (7m (S2) = 1) ||...|| (7m (stvJ = 1)) (3) 

Applying the principle of inclusion and exclusion we have: 

FK{k) = J2Phm{Si) = l)- ^((7m(s*) = 1), (7m (S,) = !)) + ... (4) 

i 

hj 
i 7^ j 

To simplify the notation, let us define W = (W (1), W (2),..., IU(|1U|)) as an ordered subset of 
(1, 2,..., Nk), where |1U| < Nk is the cardinality ofW. Now consider the sets svi/(i), sw( 2 ), sw(\w\)- 

It is clear that 

/\w\ \ 

P ((7m ('SW(I)) = l) , (7m {SW(2)) = l) , •••, (7m {^W{\W\)) = l)) = -f’ I J^Tm ^ ^ I ’ 


Let US also define: 
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IWI 

r„(wo = n^'” i^wn)) ■ (5) 

1=1 

Now, by noting that there are 2^'= possible ehoiees for W, whieh we represent by wi, W 2 , W 2 Nk, 
one ean rewrite the inelusion-exelusion expression of dH) as follows: 

2^k 

Fk (k) = (-1)'"'^'^' Pr (r^ K) = 1 ). ( 6 ) 

i=i 

Now let us look at veetor Cm defined above. There are 2^ possible realizations for Cm whieh 
are in faet the 2^ distinet binary n-tuples. We refer to these binary ?7,-tuples by Ui, U 2 ,U 2 n. 

Now the joint probability of events Tm (wj) = 1 and Cm = Ui ean be ealeulated as: 

2iv 

P{Tm (Wj) = l,Cm = Ui) = '^P {Tm (Wj) = 1, = Ui, Cm-1 = Ui) (7) 

l=l 

where we can write: 

P (Tm {^j) 1) Cm Ui, Cm—1 Ui) P (T^ (^j) 1) {Wj) 1, Cm Ui, Cm—1 Ui) . 

( 8 ) 

Note that in ([8]) 

P (Tm {u!j) 1, {Wj) 1, Cm Ui, Cm—1 Ui) 

P (Tm {u)j) 1, C Ui, Cm—1 Ui) X P {Wj) 1, jT^ (^j) 1) Cm Ui, Cm—1 Ui) 

and P (r^.i {wj) = {wj) = 1) = 1. 

Using dHl), one can obtain 

P (Tm {^j) 1) Cm Ui, Cm—1 Ui) 

P (Tm— 1 {u!j) Ij Cm—1 Ui) X P (T^ ipij) 1, Cm Ui\Tm—l {^j) I 5 Cm—1 Ui) . 

(9) 

Let us rewrite the second term in the righthand side of dH) as follows: 

P (Tm {^j) 1) Cm 1 (^j) 1) Cm—1 Ui) 

P {Cm 1 {u!j) 1, Cm—1 Ui) X P (T^ ipij) \\Cm Ui, {Wj) 1, Cm—1 Ui) . 

( 10 ) 

Note in the righthand side of dlQl) that given Cm-i, Cm is independent of T^-i {wj)- Also given 
Tm-i {wj) and Cm, {wj) is independent of Cm-i - This second claim is made by noting that if 
Tm-i {wj) = 1, then Tm {wj) = 1 if and only if given the channel realization Cm, Wj is such that 
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for every subset z = 1 : \wj\, at least one of the sensors in the subset has a Good souree- 
sensor ehannel. Therefore, it is elear that P (T^ (wj) = l\Cm = Ui, T^-i (wj) = 1, Cm-i = ui) 
is a funetion of ehannel realization n* and the set wj. If we refer to this funetion by dj (i), one 
ean write: 


dj (i) 


1 ; nls' (eL. “i (*«,(0 (ft))) > 0 

0; otherwise 


( 11 ) 


To elarify this definition, note that if Ylh=i W'j is a positive number, then given Cm = 

Ui, the subset covers the mth bit. In fact (fTTI) states that dj (i) is 1 if for every set 
at least one of the sensors in this set receives the mth bit through a Good channel. 

Let us define a matrix Q = [q{i,l)] where q{i,l) = P{Cm = Ui\Cm-i = Ui). From the 
channel model, we can observe that: 


N 

? (l 0 = n -P {Cm {n) = Ui (n) \Cm-i {n) = ui (n)) (12) 

n=\ 

and P {Cm (n) = Ui (n) |Gm-i {n) = ui {n)) is readily expressed based on the nth source-sensor 
channel state transition probabilities (/z„,e„). 

Now if we define a matrix 

~ (l 0] ) 

where Uj {i, 1) = dj {i) q {i, 1), using (fTOl) and above discussion, one can note that 

P (Tm {^j) 1) Cm 1 (^j) I 5 Cm—1 Uj (z, /) (13) 

and hence dH) can be represented as: 


P (Tm (^j) 1 ) Cm Cm—1 P (Tm—1 (^j) 1 ) Cm—1 '^l) ^ Uj (z, /) . (14) 

To simplify (fT4l) . let us define a vector 

X^ = (X^(1),X^(2),...,X^(2'^)) 


where 


(0 (^77^ I 7 C yyi • 
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Now from © and (fT4l) we have: 

2N 

Xmii) = ^aj{i,l)Xm-i{l) (15) 

1=1 

which leads to the following recursive matrix equation: 

Xm = AjXm-1- (16) 

Note that to simplify the notation, we dropped dependence of Xm to j. Also note that Aj is 
constructed by forcing some rows of matrix Q to zero. Those are the rows i such that dj(i) = 0. 
Now from (fT^ . we arrive at the following solution for Xm' 

Xm = Af-^X, (17) 

where is the m — 1 power of matrix Aj, and the initial vector Xi is expressed as: 

Xi (i) = P (Ti (wj) = 1, Cl = Ui) = P (Cl = Ui) P (Ti (wj) = l|Ci = Ui.) (18) 
It is straightforward to show that 

P {Ti{wj) = l\Ci = Ui) = dj (i) . (19) 

Now noting the independence assumption for source-sensor channels, we have: 

N 

Xi (z) = d, (z) X J] P (Cl (rz) = Ui (rz)) (20) 

n=l 

Following ifTSl we let the initial channel state Ci (rz) have the steady state probability distribution 
of the corresponding Markov process. For the Markov process of Fig. [2] this steady state 
distribution is as follows: 

P(Ci(rz) = l) = —^ (21) 

and 

P(Ci(zz) = 0) = ^^ (22) 

After solving (flSl) . we calculate F m (wj) as follows: 

2iv 

Tm {u!j) = Xm (i) 

i=l 

and by substituting in ®, we can evaluate FK{k). 


( 23 ) 
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To assess the aceuraey of our analyses, in Fig. [5] we eompare Fk (k) found using dH) with 
simulations. For these simulations, 10^ souree sequenees of length M = 128 bits are transmitted 
to = 5 sensors via identieally distributed Gilbert-Elliott channels with parameters (/i„, e„) = 
(0.0191, 0.0256) and 10® realizations of K are generated by comparing the 5 corresponding 
CSIs. The channel parameters are taken from Table U (see Section IVTl). It is clear from 

these results that our analysis is in excellent agreement with the simulated results. 

As observed from Fig. [51 Fk (5) ~ 0.5, i.e., in almost 50% of the time, employing all five 
sensors is not sufficient to cover all source bits. However, as shown in Fig.Hl in these cases, our 
algorithm forces all sensors to transmit their observations to the FC, i.e. we force Fk (N) = 1. 
In the following section, we will show that by increasing N, the coverage probability increases 
where the actual values of Fk (N) (before forcing to one) are much closer to one. 

Now, the expected value of K can be expressed as: 

N N-1 

E[K] = J2kxfK{k) = NxFK {N) -J^Fk {k) 

k=l k=l 

where by noting Fk (N) = 1 for our scheme, we reach: 

N-l 

E[K] = N-Y,Fk {k) . (24) 

A:=l 

In Section [^ we use the expected value of the subset size, E [K], to evaluate the power 
reduction achieved by our proposed algorithm. 

V. Complexity and Asymptotic Performance of the Proposed Algorithm 

In what follows, we analyze the computational complexity and asymptotic performance of the 
proposed two-phase transmission algorithm. 


A. Complexity Considerations 

Calculating Fk (k) from ® introduces a computational complexity that is exponentially 


increasing by where N^. = 


k 


Calculation of Fk (k) and consequently, E [K] is 


'Note that this computational complexity only applies to our analysis. Implementing the algorithm at the FC is considerably 
less complex as in that case the FC has the CSI realizations and only needs to compare them to find the minimum size subset. 
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time-consuming for large values of N. In fact, the run time for networks with more than N = 7 
sensors is very large. Therefore, it is desired to introduce bounds on E[K]. It is possible to 
introduce two simple upper bounds on E [K] as follows: 

Let Ek (k) be a lower bound for Ek (k). Then, from (l24l) one can find an upper bound as 
follows: 


N-l 

E[K]<N-J2^K{k). (25) 

k=l 

One possible choice for Ek (k) is by applying Bonferroni’s lower bound It20l . Let < iVfc/2 be 
an integer, then the inclusion-exclusion formula of dH) can be lower-bounded as Ek (k) > Ek (k) 
where 


2^k 

Fk (k) = (-1)'"^'^' Pr (L™ (w,) = 1). (26) 

J = 1 

\wj\ < 2Lk 

Through simulations, we concluded that for values of which introduce a reasonable compu¬ 
tational complexity, the bound of (l26l) is not tight and in fact leads to a negative value in most 
cases. 

Another simple upper bound can be derived by noting that Ek (k) > Ek (1), for k = 1 : N — 1, 
which by using (l24l) leads to: 


E{K)<N-{N-1)Ek{1) (27) 

where Ek (1) is the probability that one sensor covers the source sequence (i.e., the probability 
that at least one of the N sensors receives all M source bits via Good source-sensor channels). 
Fortunately the value of Ek (1) can be simply derived as follows. The probability that the nth 
sensor covers all source bits equals the probability that the corresponding source-sensor channel 
is initially at a Good state and stays at this state for the next M — 1 bit intervals. This probability 
is equal j (1 — e„)^~\ Therefore, the probability that none of the N sensors covers the 

source sequence equal Y[n=i and eventually, the probability that at 

least one of these N sensors covers the source sequence is given by: 
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N 


F*-(i)=i-n 


n=l 




IJ-n + 


(1 


\M-l 


(28) 


Note that if all source-sensor channels have identical parameters (/in,e„) = (/i,e), it is easy to 
verify that Fk (1) is a monotonically increasing function of N. This is expected, as by increasing 
the number of sensors, there is a higher probability that at least one of these sensors covers the 
source sequence. 

By replacing (1) from (l28l) in (1271) . we find a simple upper bound for E[K] as follows. 


< iV - (iV - 1) 1^1 - J] (^1 - (1 - j . (29) 

Figure [6] shows E [K] as a function of N for a network with identical source-sensor channel 
parameters (/i, e) = (0.0041,0.0075) and source sequence lengths of M = 200, 256, 300 bits. 
The values of (/r, e) are based on Table H Note the non-monotonic behaviour that is observed in 
Fig. [^for the upper bound of E[K]. This non-monotonic behaviour is due to two reasons. The 
first is based on the fact that this upper bound is not tight. However, there is another rational 
behind the non-monotonic behaviour of this upper bound. That is, in cases where the source is 
not covered by any subset, we demand transmission from all N sensors (i.e., K = N). When N 
increases, there are subsets with larger sizes to examine for possible coverage. Therefore, E [K] 
might increase in such cases. Although the upper bound of (l29l) is not tight, as we will see in 
Section |Vll even by applying this simple bound, we observe considerable power reduction when 
employing our proposed algorithm for networks with large values of N. 


B. Asymptotic Performance 

Here, we consider the asymptotic performance of our proposed algorithm for large values of 
N. For simplicity, let us assume that all N source-sensor channels have identical state transition 
probabilities (/i,e). It is clear from (|2^ that for identical values of (/in,en) = (/i,e), Fk (1) is 
a monotonically increasing function of N and lim^^oo (1) = 1. By noting that TV (1) < 
Fr (/c) < 1 for fc = 2 : N, the lower bounds of Fr {k) = Fr {!), k = 2 : N dxo asymptotically 
tight. Therefore, the upper bound of (l29l) is asymptotically tight. If we let (/i„, e^) = {p, e) and 
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by taking the derivative of (l29l) with respeet to N, one can find a value Nq such that for all 
N > Nq, the upper bound of E[K] is monotonically decreasing by N 


No = 


1 + 



(30) 


where x = (1 “ From the above discussion, we conclude that for sufficiently 

large values of N, is a monotonically decreasing function of N (decaying by a rate of 
or faster). Rewrite O as 


1 ^ P I E[K] 
T] - N N 


(31) 


and note that p < 1. We observe that 3, which is the ratio of power consumption for our 
proposed algorithm to the conventional transmission scheme, decays by increasing N (at least 
by a rate of ;^). Therefore, our proposed algorithm becomes asymptotically more power efficient 
by increasing N. 

At the end of this section, we note that when we were motivating the idea in Section H we 
applied a parameter u for bounding the distortion, where we defined u as the probability that 
there exists a (minimum-size) subset that covers all source bits. We claimed that for sufficiently 
large N, v can be arbitrarily close to one. To prove this claim, note that the probability that such 
subset exists, is greater than or equal the probability that such subset exists and its size is k (for 
an arbitrarily chosen k < N). Therefore, u > Fk (k) > Fk (1) and F^ (1) could be arbitrarily 
close to one, given a sufficiently large N. 


VI. Numerical Results 

In this section, we provide some numerical results to evaluate the power efficiency of our 
proposed algorithm. In this work, and without loss of generality, we only consider cases where all 
source-sensor channels have identical parameters (p,e). The Gilbert-Elliott channel parameters 
are derived by simulating a correlated Rayleigh fading channel using Jakes model and then 
quantizing the simulated channel by assuming a threshold on the fading amplitude. If we represent 
the fading amplitude by a and assume that the source is transmitting each bit with energy Ff, 
and the AWGN has a one-sided power spectral density of A"o, then the instantaneous received 
SNR equalsat the sensor. We consider a threshold of athr = 1- That is we assume that 
SNRs above ^ leads to a Good delivery of the source bit to the sensor (i.e, the probability of 
detection error, pc, is sufficiently low to have pc < D as discussed in Section U). The assumption 
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Fig. 5. The cumulative mass function of selected subset size, K, for a network with N = 5 sensors and source sequence length 
of M = 128 bits. The source-sensor channels have identical transition probabilities of (/in, £n) = (0.0191, 0.0256). 



Fig. 6. The upper bound on E[K] for a network with identical source-sensor channel parameters (/i, e) = (0.0041,0.0075). 
The source sequence lengths are M — 200, M = 256, and M = 300 bits. 
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of athr = 1 is justified as follows. If we assume that the channel phase shift is perfectly estimated 
and compensated at the sensor node, then for all a > athr, the channel provides error detection 
probabilities less than or equal to the error detection probability of an AWGN channel with SNR 
of Therefore, by setting this threshold, we eliminate the non-constructive effect of fading 
and provided source-sensor channels with link qualities equivalent or superior to an AWGN 
channel. We consider a slow-fading channel, i.e., channels with the normalized fading rates of 
fdTs < 0.01, where fd is the maximum Doppler shift and is the symbol duration. The reason 
we consider slow-fading channels is that as discussed in previous sections, for these channels 
the run-length coding of CSI sequences provides an efficient compression. 

As discussed earlier, the parameters (/r, e) for the Gilbert-Elliott channel are estimated using 
Monte-Carlo simulation of sufficiently large number of realizations of the fading channel am¬ 
plitude. For values of fdTs = 0.002, 0.005, 0.008, the corresponding values of (/i, e) are shown 
in Table H Table H] also shows the expected value of compression rate for CSIs, for different 
sequence lengths of M = 128, M = 256. As expected, the compression rate decreases when 
increasing M. 

Tables H] and In] show values of {E [K] , r/) for networks with N = 4, 5, 6 sensors. We observe 
that E [K] is a non-monotonic function of N. The justification of this non-monotonic behaviour 
was discussed in Section IV-BI and Fig. 0 Note that the efficiency factor, rj, is monotonically 
increasing function of N, which confirms the increase in efficiency of our proposed algorithm 
as the number of sensors, N, increases. 

From Tables jn] and jnll it is clear that our algorithm is more efficient for channels with slower 
fading rates. For instance, in Table HI if we let N = 5, we observe that the efficiency factor 
for channels with fdTg = 0.002 is 1.92 which shows an almost two-fold decrease in power 
consumption achieved by our algorithm compared to the conventional transmission scheme. 
However, when we increase fdTs to 0.008, r] decreases to 1.30. Also, by comparing the results 
of Table jn] and Table HIl we observe that our proposed algorithm is more efficient for shorter 
source sequence lengths, M. The reason is that we defined a coverage event as the event that all 
source bits are covered. Therefore, the coverage probability of a subset decreases by increasing 
M. This results in an increase of E [K] and consequently a decrease in rj. The only case where our 
scheme shows an inferior performance to the conventional transmission scheme is for M = 256, 
fdTs = 0.008, and = 4, where rj = 0.99. 

To examine the case of large networks with large values of N, we turn to the upper bound 
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TABLE I 

Gilbert-Elliott channel transition probabilities and achieved compression rates by run-length coding 

SCHEME FOR DIFFERENT VALUES OF THE NORMALIZED FADING RATE, fdTs. THE FADING AMPLITUDE THRESHOLD FOR 
DECIDING BETWEEN GOOD AND BAD STATES IS SET TO 1. THE SOURCE SEQUENCE LENGTHS OF M = 128 AND M = 256 

BITS ARE CONSIDERED. 


fdTs 

e 


p{M = 128) 

p{M = 256) 

0.002 

0.0075 

0.0041 

0.1071 

0.0813 

0.005 

0.0165 

0.0112 

0.1630 

0.1454 

0.008 

0.0256 

0.0191 

0.2223 

0.2134 


TABLE II 

Values of(E [K] , jy) for networks with N = 4,5,6 sensors and source sequence length ofM = 128 bits. 


fdTs 

N = 4 

iV = 5 

iV = 6 

0.002 

(2.44,1.57) 

(2.49,1.92) 

(2.43, 2.37) 

0.005 

(2.97,1.28) 

(3.12,1.52) 

(3.05,1.87) 

0.008 

(3.33,1.13) 

(3.63,1.30) 

(3.55,1.59) 


of (l29l) . Replacing this upper bound in ([I]) provides a lower bound on rj. Figure |7] shows the 
upper bound of E[K] and the resulting lower bound on r] for networks with source sequence 
length of M = 256 bits and f^Tg = 0.002. One can note the considerable gains for these large 
values of N when using our algorithm. For example, for a network with iV = 50 sensors, our 
proposed algorithm provides at least a twelve-fold decrease in the consumed power compared 
to the conventional transmission scheme with all nodes transmitting {t] > 12). To examine the 
effect of different block lengths on rj, we also consider block lengths M = 200 and M = 300 
in Fig. Ill As observed, the efficiency factor decreases by increasing the block length. This is 
due the fact that as M increases, the probability that k sensors cover all M bits decreases. As a 
result E [K] increases and rj becomes smaller. Nonetheless, we observe that for M = 300 and 
N = 50, our algorithm has an efficiency factor close to 6. 

VII. Conclusion 

We analyzed a WSN where source-sensor channels are modeled as quantized correlated fading 
channels (Gilbert-Elliott channels). We proposed a two-phase transmission scheme where at the 
first phase compressed channel state information sequences are transmitted to the FC and a 
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TABLE III 

Values of{E [K] , rj) for networks with N = 4,5,6 sensors and source sequence length oeM = 256 bits. 


fdTs 

N = 4 

N = 5 

N = 6 

0.002 

(3.06,1.27) 

(3.25,1.50) 

(3.31,1.77) 

0.005 

(3.62,1.06) 

(4.05,1.19) 

(4.23,1.37) 

0.008 

(3.83,0.99) 

(4.46,1.07) 

(4.78,1.20) 



Fig. 7. Values of rj for networks with source sequence lengths ofM = 200, 256, 300 bits and fdTs = 0.002. The dashed line 
shows E \K] for M = 256 bits. 


subset of sensors are seleeted to transmit their observations to the FC at the seeond phase. Also, 
we analytieally derived the probability distribution of the size of the seleeted subset and the 
expeeted value of this subset size. We presented simulation results to assess the accuraey of our 
analyses. We defined an effieieney faetor for our proposed algorithm and evaluated this faetor for 
several ehannel eonditions and network setups. In most oases our proposed two-phase algorithm 
showed a superior power effieieney compared to a conventional one-phase transmission scheme 
over slow-fading channels. 
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